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Abstract 

We consider the two-dimensional randomly site diluted Ising model and the random-bond 
±J Ising model (also called Edwards- Anderson model), and study their critical behavior at the 
paramagnetic-ferromagnetic transition. The critical behavior of thermodynamic quantities can be 
derived from a set of renormalization-group equations, in which disorder is a marginally irrele- 
vant perturbation at the two-dimensional Ising fixed point. We discuss their solutions, focusing in 
particular on the universality of the logarithmic corrections arising from the presence of disorder. 
Then, we present a finite-size scaling analysis of high-statistics Monte Carlo simulations. The nu- 
merical results confirm the renormalization-group predictions, and in particular the universality of 
the logarithmic corrections to the Ising behavior due to quenched dilution. 

PACS numbers: 75.10.Nr, 64.60.Fr, 75.40. Cx, 75.40.Mg 
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I. INTRODUCTION AND SUMMARY 



Random Ising systems represent an interesting theoretical laboratory in which one can 
study general features of disordered systems. Among them, the two-dimensional (2D) 
random-site and random-bond Ising models have attracted much interest. In particular, 
the effects of quenched disorder on the critical behavior at the paramagnetic-ferromagnetic 
transitions, which are observed for sufficiently small disorder, have been much investigated 
and debated, see Refs. 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27, 
28,29,30,31,32,33,34,35,36,37,38,39. 

Renormalization-group (RG) and conformal field theory 2,4 ' 8 ' 12 predict the marginal ir- 
relevance of random dilution at the paramagnetic-ferromagnetic transition. Therefore, the 
asymptotic behavior is controlled by the standard Ising fixed point, characterized by the 
critical exponents v — 1 and r\ = 1/4; disorder gives only rise to logarithmic corrections. 
The marginality of quenched disorder coupled to the energy density, as it is the case for ran- 
dom dilution, is already suggested by the Harris criterion, 40 which states that the relevance 
or irrelevance of quenched dilution depends on the sign of the specific-heat exponent of the 
pure system; in the case of the 2D Ising model, the specific heat diverges only logarithmically 
at the transition, i.e. a = + . The marginal irrelevance of disorder has also been supported 
by numerical studies of lattice models; see, e.g., Refs. 10,23,24,26,28,32,33,34,36,37,38 (see, 
however, Refs. 13,14,30 for different scenarios). We recall that in three dimensions ran- 
dom dilution is a relevant perturbation of the pure Ising fixed point, leading to a new 
three-dimensional randomly diluted Ising (RDI) universality class, which is characterized by 
different critical exponents; see, e.g., Refs. 41,42. 

In this paper we return to the issue of the critical behavior of 2D randomly diluted Ising 
systems. By using the RG results reported in Refs. 3,5,6,28, we show that their critical 
behavior can be derived from the RG equations 
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where / is the flow parameter (the logarithm of a length scale), uj, u t , and Uh are the scaling 
fields associated with the leading operators of the three different conformal families of the 
2D Ising model, i.e. the identity, energy, and spin families, and g is the marginally irrelevant 
scaling field associated with disorder. Higher-order terms in Eqs. (1) are not necessary, 
because they can be reabsorbed by appropriate analytic redefinitions of the scaling fields. 
The appearance of the term diu 2 in the first equation, where dj is a nonuniversal constant, 
is due to the resonance of the identity operator with the thermal operator, which already 
occurs in the pure Ising model. 43 It is interesting to note that randomness couples only to 
the thermal scaling field u t . It would be interesting to understand if these conclusions also 
apply to the irrelevant operators, i.e., if the only operators that couple disorder are those 
that belong to the conformal family of the energy. 

The analysis of the RG equations shows that random dilution gives rise to logarithmic 
corrections which are universal after an appropriate normalization of the scaling field as- 
sociated with disorder. Additional scaling corrections due to the irrelevant operators are 
suppressed by power laws as in standard continuous transitions. For these reasons, we pre- 
fer to distinguish the randomly dilute Ising (RDI) critical behavior characterized by the RG 
equations (1) from the standard 2D Ising universality class of pure systems, although they 
share the same 2D Ising fixed point. 

The RG equations (1) allow us to determine the scaling behavior of any thermodynamic 
quantity. Denoting with t, h, p, and L the reduced temperature, the magnetic external field, 
the disorder parameter, and the lattice size, respectively, the free energy satisfies the scaling 
equation 



(we consider models defined on square L x L lattices with periodic boundary conditions), 
where «/(/), u t (l), u h (l), and g{l) are the solutions of the RG equations. From Eq. (2) one 
can derive the scaling behavior of the relevant thermodynamic quantities and determine the 
logarithmic corrections due to the quenched disorder. At the critical point t — h — 0, we 
obtain the asymptotic behaviors 28 



F(t, h, p, L) = e- 2l uj(l) + e- 2l f(u t (l), u h (l), g(l), e l L) 



(2) 



Ch ~ In In L 



(3) 



for the specific heat, and 
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/?= R"Mu(hiL))= R fO(j^). (5) 
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^ ' \/Id.L 
for the magnetic susceptibility x> an y RG invariant quantity i?, such as the quartic Binder 
cumulant £/ 4 and the ratio = £/L, and its derivative with respect to the temperature. 
Here R* indicates the Ising fixed-point value and <?(mL) is the solution of Eq. (1) with 
I = InL. For L — > oo, <?(lnL) behaves as 



gQnL) ' 



In lnL/Ln 

1 H — - + 
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where L is a length scale. The functions f x (x), fii(x), and fdR^x) are normalized such that 
/#(0) = 1 and are universal once the scaling field g(\nL) is appropriately normalized. In 
the above-reported equations we have neglected scaling corrections which are suppressed by 
power laws. They are due to the analytic dependence of the scaling fields on the Hamilto- 
nian parameters, to the background (i.e., the contribution of the identity operator in the 
RG language), and to the irrelevant operators. 44 ' 45 In particular, we expect scaling correc- 
tions associated with the leading irrelevant operator appearing in the pure Ising model (the 
corresponding exponent is oo = 2). 

Moreover, in this paper we compare the theoretical predictions with a finite-size scaling 
(FSS) analysis of numerical Monte Carlo (MC) results for the randomly site-diluted Ising 
model and for the random-bond ±J Ising model, also known as Edwards- Anderson model. 
Our main results can be summarized as follows. Our FSS analyses provide a robust ev- 
idence that the paramagnetic-ferromagnetic transitions in these models present the same 
RDI critical behavior. Note that this implies that frustration in the random-bond ±J Ising 
model is irrelevant at the paramagnetic-ferromagnetic transition line. The FSS behaviors 
are in agreement with the predictions of the RG equations (1). The asymptotic critical 
behavior appears to be controlled by the Ising fixed point. The logarithmic corrections and 
their universal behavior are consistent with the theoretical results obtained from Eqs. (1), 
cf. Eqs. (2), (3), (4), (6), (7). 

The paper is organized as follows. In Sec. II we define the randomly site-diluted model 
and the ±J Ising model, we briefly discuss their phase diagrams, and define the quantities 
that we consider in the paper. In Sec. Ill we discuss the RG flow at a 2D Ising fixed point in 
the presence of a marginally irrelevant operator, and the implications for the infinite- volume 



and finite-size critical behavior. In particular, we focus on the universal features of the 
logarithmic corrections due to disorder. Finally, in Sec. IV we present our FSS analysis of 
high-statistics MC results for the randomly site-diluted and random-bond ± J Ising models. 

II. RANDOMLY SITE-DILUTED AND RANDOM-BOND ±J ISING MODELS 
A. The models and their phase diagrams 

The randomly site-diluted Ising model (RSIM) is defined by the Hamiltonian 

H s = - ^2 P x Py a x°y> (8) 

<xy> 

where the sum is extended over all pairs of nearest-neighbor sites of a square lattice, a x — ±1 
are Ising spin variables, and p x are uncorrelated quenched random variables, which are equal 
to one with probability p (the spin concentration) and zero with probability 1 — p (the 
impurity concentration). The RSIM is expected to undergo a continuous transition for any 
V > Ppcro where 46 p perc = 0.59374621(13) corresponds to the site-percolation point of the 
impurities; moreover, T c — > for p — > p pe rc, see, e.g., Ref. 47. For p < p perc the ferromagnetic 
phase is absent. Thus, the paramagnetic-ferromagnetic transition line starts from the pure 
Ising point X Is = (p = 1,T = T ls ), where T ls = 2/ln(l + y/2) = 2.26919... is the critical 
temperature of the 2D Ising model, and ends at X pcrc — {p — p pcTC ,T = 0). Along this 
line the critical behavior is expected to be universal, i.e. independent of dilution, and to be 
characterized by the RG equations (1). As we shall see, this is supported by the analysis of 
our MC results. 

The random-bond ± J Ising model, also known as Edwards-Anderson model, 48 is defined 
by the lattice Hamiltonian 

Hb = - ^ Jxy&x&y, (9) 

(xy) 

where a x = ±1, the sum is over all pairs of nearest-neighbor sites of a square lattice, and 
the exchange interactions J xy are uncorrelated quenched random variables, taking values 
±J with probability distribution 

P(Jxy) = pt(Jxy ~ J) + (1 " p)S{J xy + J). (10) 

In the following we set J = 1 without loss of generality. For p = 1 we recover the standard 
Ising model, while for p — 1/2 we obtain the bimodal Ising spin-glass model. The ±J Ising 
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FIG. 1: Phase diagram of the square-lattice random-bond ±J Ising (Edwards- Anderson) model 
in the T-p plane. 

model is a simplified model 48 for disordered spin systems showing glassy behavior in some 
region of their phase diagram. Its phase diagram in two dimensions is sketched in Fig. 1 (it 
is symmetric for p —>■ 1 — p) . For sufficiently small values of the probability p a = I — p of 
the antiferromagnetic bonds, the model presents a paramagnetic phase and a ferromagnetic 
phase. The paramagnetic-ferromagnetic transition line starts from the Ising point Xj s = 
(p = 1,T = Tj s ) and extends up to the multicritical Nishimori point (MNP) X* = (p*,T*), 
located along the so-called Nishimori line (N line) defined by 2p — 1 = Tanh(l/T), 49 ' 50 ' 51 ' 52 
with 53 p* = 0.89081(7) and T* = 0.9528(4). The critical behavior is expected to be in 
the same universality class as that of the transition in the RSIM. As we shall see, our FSS 
analysis strongly supports this scenario. A detailed discussion of the phase diagram can be 
found in Ref. 53 and references therein. 



B. Observables 

In our FSS analyses we consider models defined on a square L x L lattice with periodic 
boundary conditions. The two-point correlation function is defined as 

G(x) = [(a a x )], (11) 
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where the angular and square brackets indicate the thermal average and the quenched aver- 
age over disorder, i.e. over p x in the case of RSIM and over J xy in the case of the ± J Ising 
model. We define the magnetic susceptibility x = XL G(x) and the correlation length £, 

e ^ G(0)-G(g min ) ; (12) 

?min^'(?min) 

where g min = (2tt/L,0), q = 2 sing/2, and is the Fourier transform of G(x). We also 
consider quantities that are invariant under RG transformations in the critical limit. Beside 
the ratio 

% = tlL, (13) 
we consider the quartic cumulants U4 and U 2 2 defined by 

Ua = t — 12' U22 = — f — 12 — ' 

where 

V k = ((J2^) k )- (15) 

X 

The above RG invariant quantities R^, U4, and U22 are also called phenomeno logical cou- 
plings. In the critical (T = T c ) 2D pure Ising model, they converge for large L to the 
universal values 54 

R* = R ls = 0.9050488292(4), (16) 

U* = Ui s = 1.167923(5), (17) 

U* 2 = 0. (18) 

Finally, we consider the derivatives 

ff ^f> * = (19) 

which can be computed by measuring appropriate expectation values at fixed (5 and p. 



III. RENORMALIZATION-GROUP FLOW AND FINITE-SIZE SCALING 

A. Ising RG flow in the presence of a marginally irrelevant scaling field associated 
with disorder 

In this section we discuss the RG flow close to the 2D Ising fixed point in the presence 
of a marginally irrelevant scaling field associated with disorder. 
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Let us consider a system with a marginal scaling field u = g and with a set of scaling 
fields iik, k > 1, with RG dimensions y k 7^ 0. Close to the fixed point # = u\ — . . . — 0, the 
RG equations have the form 

^ = ^ \ijUiUj + ^ b 0,ijmUiUjU m + ■■■ 
0<i<j 0<i<j<m 

= VkUk + ^ b k,ijUiUj + ^2 KijmUiUjUm + ... (20) 
0<i<j 0<i<j<m 

If there are no degeneracies (yk ^ Uh for all k ^ h) and no resonancies (i.e., there is no 
combination with integer coefficients of the RG dimensions that vanishes), one can redefine 
the scaling fields in such a way to simplify the RG equations. We define 

9 = 9 + ^2 C 0,^% + C 0,ijmUiUjU m + ... (21) 

0<j<j 0<i<j<m 

U k =Uk+ ^2 C k,ijUiUj + ^2 C k,ijmUiUjU m + ... (22) 
0<i<j 0<i<j<m 

With a proper choice of the coefficients c^-, Ck,ij m , ■ ■ ., we can simplify the RG equations, 
obtaining the simple canonical form: 

| = -M 2 - hg\ (23) 
—jj- = VkUk + c k gu k . (24) 

By normalizing appropriately the scaling field g, we can also set 1 62 1 — 1 ■ in the case we 
are considering, g is marginally irrelevant so that b 2 > (we assume that g is defined such 
that g(l — 0) > 0). We can thus set b 2 = 1. Once this choice has been made, 63 and all 
coefficients c k are universal. 

The simple form we have derived above does not strictly apply to the RSIM. Indeed, in 
the Ising model the RG operators belong to three different conformal families and within 
each family the RG dimensions differ by integers (see Ref. 45 for a discussion of the irrelevant 
operators in the pure Ising model). Thus, in the present case there are both degeneracies 
and resonancies. If we limit our considerations to the relevant scaling fields, we must only 
consider the resonance between the identity operator and the thermal operator, which is 
responsible for the logarithmic divergence of the specific heat in the pure Ising model. 43 By 
a proper redefinition of the nonlinear scaling fields, one can show that in this case the RG 
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equations (20) for the relevant scaling fields can be written as: 



duj 



2uj + cjgui + diuf, 



(25) 



dl 

du t 

~dT 

du h 



VhU h + c h gu h , 



y t u t + c t gu t , 



(27) 



(26) 



dl 
dg 

dl 



-g 2 - hg\ 



(28) 



where the couplings to the irrelevant scaling fields due to the additional resonancies have 
been neglected. The scaling field uj is associated with the identity operator. The additional 
term dju^ which appears in Eq. (25) is due to the resonance with the thermal operator, 
as in the pure Ising model. 43 The scaling fields Ut and Uh are the relevant scaling fields 
associated with the reduced temperature t and the external field h, respectively; y t — 1 
and Uh— 15/8 are the corresponding RG dimensions. Finally, g is the marginally irrelevant 
operator associated with randomness. The coefficients cj, c t , Ch, and 63 are universal, being 
independent of the normalization of the scaling fields. By using conformal field theory, c t , 
Ch, and 63 have been computed: 3 ' 5 ' 6,19 ' 28 



Let us now integrate the RG equations. Since 63 < 0, Eq. (28) has two fixed points with 
g > 0: one is g = and is stable; the second one is g — — l/fo 3 = 2 and is unstable. Thus, 
the basin of attraction of the Ising FP corresponds to g = g(l = 0) < — 1/63 = 2; for 
go > 2, g(l) flows to infinity. It is important to note that Eq. (28) is only valid within the 
basin of attraction of the stable fixed point g = 0. The redefinitions of the scaling fields 
that we have used to obtain the simple canonical form (28) cannot be extended outside the 
basin of attraction since they are expected to become singular at the unstable fixed point. 
The presence of an unstable fixed point indicates that the behavior for large values of the 
disorder is not controlled by the Ising fixed point. The RG flow could be attracted by a new 
fixed point — thus, for large values of the disorder the transition would be continuous and 
in a new universality class — or could go to infinity, indicating the absence of a continuous 
transition. A similar phenomenon was conjectured in three dimensions 55 on the basis of a 
perturbative field-theoretical analysis of the RG flow. 

If go < — 1/&3, the function g(l) is implicitly given by (we do not replace 63 with its 



<h = -1/2, 



c h = 0, 



1/2. 



(29) 
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theoretical value —1/2, in order to obtain general expressions that can be tested numerically) 



F[g(l)]=l + F(g ), 

F(x) = - + b 3 ln(—— 

x \l + b 3 x J 



(30) 



The solution can be simplified if we introduce 



9(1) 



i + MO' 



(31) 



which satisfies the implicit equation 



FW)]=l + F(9o), 



Fix) = — h b 3 Inx. 

x 



(32) 



This equation can be inverted to give 



(33) 



The function $(x, I) cannot be computed analytically. However, it is easy to determine it 
in the large-/ limit. We obtain 



1 b 3 kiy b 2 3 (ln 2 y - lay) 



+ 



b\ In 3 y 



y r r v yr 

where y = I + F(g ). Since Ch = the equation for Uh gives 

u h (l) = u hfi e Vhl , 



(34) 



(35) 



where Uh t o — u^{l — 0). In order to determine u t (l), we rewrite the corresponding equation 

as 

du t gdg dg 

= y t dl - (h = y t dl - c t — , (36) 

ut 9 2 + b 3 g 4 g 



which gives (q = -1/2, y t = 1) 



u t {l) = u tfl e l 



m 

9o 



1/2 



(37) 



where u t:0 — u t (l — 0). For large / the function u t (l) behaves as 



u t (l) = u tt0 g l,2 ^j2 



(38) 
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Let us finally consider the identity operator. If dj 
the case of u t . Thus, we write 



the solution can be obtained as in 



ui(l) = u Ifi e 2l K(l) 



m 

9o 



1 -ci 



where K{1) is an unknown function of /, such that K{1 = 0) 
equation for Ui(l) and using the result for u t {l), we obtain 



dK 
~dT 



d r 



u 



Ul,0 



and therefore 



ffl 



cj+1 



K{1) = 1 - d^-g^'- 1 / dxx c '-\l - b 3 x). 
ui, Jg 



(39) 

Substituting in the 

(40) 
(41) 



The behavior of uiil) for / — > oo depends on the value of c/. Since the integral appearing in 
K(l) diverges as l~ Cl for cj < 0, as In / for c/ = 0, and is finite for c/ > 0, we obtain 

{e 2l l c ' for a > 0, 

e 2l \nl for C/ = 0, (42) 
e 21 for c 7 < 0. 

The RG equations do not fix completely the normalization of the scaling fields. First, 
one can redefine u t , Uh, and ui by a multiplicative constant; 56 such a redefinition is not 
possible for g, since a multiplicative constant would break the condition b 2 = 1. Beside 
these trivial redefinitions there is also a nonlinear set of transformations that leave the 
equations invariant. Given a constant A, we define g\ as the solution of the equation 



Then, for any A we have 



F(g x ) = F(g) + A. 



d 9\ 2,3 



(43) 



(44) 



Note that the transformation is analytic in a neighborhood of g — 0. If g\ is defined as in 
Eq. (31), we obtain 

~g x = g[l-X~g + 0(g 2 )}. (45) 

Analogously, if we define 

u t ,x = u t (~g/~g x r 1/2 , (46) 

then u tt \ satisfies the same equation of u t with g x replacing g, as it can be seen from Eq. (37). 
A similar redefinition can be made for ui. This invariance implies that, beside fixing the 
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normalizations of u t , Uh, and ui, we must also appropriately fix g. In practical terms, F(g ) 
is completely arbitrary and must be fixed in order to define g{l) unambigously. Finally, note 
that there are no analytic redefinitions of g that map Eq. (28) in an identical equation with 
b' 3 7^ b 3 , proving the universality of b 3 . 

Neglecting scaling corrections that are suppressed by power laws, we write the free energy 
in the scaling form 43 



H^Kp) = e- 2l Ul (l)+e- 2l f sing (u t (l),u h (l),g(l)) : 



21 



(47) 



for any /. Note that the whole dependence on t e T/T c — 1, h, and p is encoded in the 
constants go, u tj o, Uh,o, ui t o, and di. Of course, u tt o ~ t and Uh,o ~ h vanish at the critical 
point, while g Q vanishes for p = 1. The independence of Eq. (47) on I allows us to simplify 
the general expression for the free energy. We choose I such that u t (l) = 1 and thus 

'In I In r I 



= r-^-lnr) 



1/2 



1 + 



lnr 



m = 



lnr 



1 + 



In | lnr | 
lnr 



(48) 



— 1/2 

where r = u t! o/g Q ' . Substituting these expressions in Eq. (47), we obtain the general 
dependence of the free energy on t and h. 

In order to determine cj, we consider the specific heat Ch- The leading singular behavior 
is due to the temperature dependence of the scaling field ui. Using Eq. (42) we obtain 

(In l/t) Cl for a > 0, 



C h oc 



dt 2 



lnln(l/t) 
constant 



for ci = 0, 
for cj < 0. 



(49) 



The asymptotic behavior of the specific heat of 2D randomly diluted Ising systems has been 
determined in Refs. 1,2,4, obtaining 



0,~lnln(l/t). 

Comparing with Eq. (49), we obtain cj = 0. In this case we have 



ui(l) 



ur e 21 - 



go 



i 9(1) , ( ~ ~ x 
In— b 3 {g- g ) 

go 



(50) 



(51) 



It is interesting to note that, since Ch = Cj = 0, randomness couples only to the thermal 
scaling field u t . This result appears quite natural from the point of view of the Landau- 
Ginzburg- Wilson approach to critical phenomena. Indeed, in field theory randomly diluted 
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models are obtained by coupling disorder to the energy operator: 57 ' 58 

H = Jd d x ^-{d^{x)f + ^(x) 2 + ^(x)<P(x) 2 + ±g [0(:r) 2 ] 2 } , (52) 

where r oc T — T c , and ip(x) is a spatially uncorrelated random field with Gaussian distribu- 
tion. The 2D RDI critical behavior has been also investigated by using this field-theoretical 
approach and the so-called replica trick. The analysis of the corresponding five-loop per- 
turbative expansions 7 ' 34 gives results which are substantially consistent with the marginal 
irrelevance of disorder. 



B. Finite-size scaling 

Let us now discuss the implications of the above RG analysis for the FSS of thermo- 
dynamic quantities at the critical point. We start from the scaling behavior of the free 
energy 

F(t, h, p, L) = e- 2l uj(l) + e- 2l f(u t (l),u h (l),g(l), e l L^), (53) 

where the contributions of the irrelevant scaling fields have been neglected. By choosing 
/ = InL, we can write 

F(t,h,p,L) = L- 2 Ul (\nL) + L- 2 f(u t (\nL),u h (\nL),g(\nL)). (54) 

If we set F(g ) = — lnL in Eq. (32), for L — > oo we have 

In In L/ Lq 



Sf(lnL) = 



In L/L, 



o 



1 " h : , / " + O 



lnL/L \(\nL/L ) 



(55) 



The free energy can then be written as 



^(tjhjPjL) = ki lnp(lnL) + k 2 + k 3 g(\nL) 

+f(u t:0 Lg(\nL) 1 / 2 ,u h:0 L 15 / 8 ,g(\nL)). (56) 

The constants k iy u tj0 , W/^o, and L Q depend on t, h, and p. Moreover, u t>0 ~ t and u h $ ~ 
u h ~ h close to the critical point. The terms proportional to ki, k 2 , and k 3 are due to the 
identity operator, whose dependence on g(lnL) is given in Eq. (51). Eq. (56) is valid up to 
contributions of the irrelevant operators, which are expected to scale as inverse powers of L. 



13 



From Eq. (56) we can compute zero-momentum quantities that involve disorder averages 
of a single thermal average. For instance, for the specific heat at T = T c and h = we 
obtain 

Cft-lnlnL. (57) 
For the susceptibility at h = we obtain 

X= (^f J L^f x (u tfi L~g(\nLy/*,~g(\nL)), (58) 

where, as before, we neglect power-law scaling corrections. A similar scaling equation holds 
for U 4 : 

U A = /^> t , Ls(lnL)V 2 ,£(lnL)). (59) 

The determination of the scaling behavior of U22 and = £/L requires an extension of the 
scaling Ansatz (56). A detailed discussion is presented in Sec. 3.1 of Ref. 42. It shows that 
both quantities behave as C/4, apart from scaling corrections. Thus, if R = U22 or R — R^, 
we also have 

R = f R (u tfi L~g(\nL) l l 2 ,~g(\nL)). (60) 

Derivatives of the phenomenological couplings have a simple behavior as well, the leading 
term being of the form 

W = (^) L ^ lnL ) V2 ^Ko^(lnL) 1/2 ^(lnL)). (61) 
At the critical point we can set u t>0 = 0, so that we can write the scaling behaviors 

R = g R [g(\nL)], 

^ = (%^)^ (lnL)1/2 ^ [ ^ (lnL)] - (62) 

The functions gn(x) and gdR^x) are universal once an appropriate normalization is chosen 
for ^(lnL), which is independent of the model. For this purpose, let us consider a phe- 
nomenological coupling R. For L — > 00 we can expand 

R = R* + ng(\n L) + r 2 ^(ln L) 2 + • • • (63) 

Now we normalize ^(lnL) by requiring r 2 = 0. It is easy to prove that this is a correct 
normalization condition. Indeed, imagine that g(lnL) has been normalized arbitrarily so 
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that r 2 7^ 0. Then, redefine g(\nL) by using Eq. (45). By properly choosing A, it is easy to 
see that one can set r 2 = 0. This condition fixes uniquely the scale L . 

Note that, in the pure Ising model, we have U 22 = 0, so that we expect at the critical 
point 

U 22 ~~g{\nL) (64) 

for L — > oo. It is natural to invert this relation to express g(lnL) in terms of U^L). Then, 
we obtain the scaling forms 



R(L) = f R (U 22 ), (65) 
X (L) = d x L 7 ^f x (U 22 ), (66) 
d dR LXj)l 2 2 ] dR (XJ 22 ), (67) 



dp 

where f R (x), f x ( x ), fdR.(x) are universal scaling functions that are normalized such that 
/r(0) = R*, /x(0) = /rfft(O) — 1 an d have a regular expansion in powers of x. Note that 
these scaling equations are much simpler than those in terms of g(lnL), since they are 
independent of the scale L and of the normalization of g(\nL). 



C. Finite-size scaling at a fixed phenomenological coupling 

Instead of computing the various quantities at fixed Hamiltonian parameters, one may 
study FSS keeping a phenomenological coupling R fixed at a given value Rf, as proposed in 
Ref. 59 and also discussed in Refs. 42,60. This means that, for each L, one considers f3f(L) 
such that 

R(J3 = P f (L),L) = R f . (68) 

All interesting thermodynamic quantities are then computed at f3 — /3f(L). The pseudocrit- 
ical temperature (3f(L) converges to f3 c as L — > oo. 

In the next section we report a FSS analysis of MC simulations keeping the phenomeno- 
logical coupling R^ fixed. The value Rf can be specified at will as long as it is between 
the corresponding high-temperature and low-temperature values. Since we wish to check 
the hypothesis that the asymptotic critical behavior is governed by the Ising fixed point, 
we choose R^j = Ri s , where Ri s = 0.9050488292(4) is the universal value of R^ = £/L at 
the critical point in the 2D Ising universality class 54 for square L x L lattices with periodic 
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boundary conditions. Note, however, that this choice does not bias our analysis in favor of 
the Ising nature of the transition. By fixing to the critical Ising value, we can perform 
the following consistency check: if the transition belongs to the Ising universality class, then 
any other RG-invariant quantity must converge to its critical-point value in the Ising model. 
In the (t, L) plane, the line = Ri s is obtained by solving the equation 

f R ^u tfi L~g(\nL) l l 2 ,~g(\nL)) = R ls . (69) 

It gives a relation 

u t ,oL£(lnL)V 2 = £;(5(lnL)), (70) 

where k(x) has a regular expansion in powers of x. Moreover, since we have chosen R^j = 
Rj s , we have k(0) = 0. Substituting relation (70) in the above-reported scaling equations 
for the susceptibility x, the phenomenological couplings R, and their derivative, we obtain 
at fixed i?£ 

x(L) = c x L^C x (g(^L)), (71) 
R(L) = C R (~g(\nL)), (72) 
= c dR L~g{\nLy/ 2 C dR {g{\nL)). (73) 

The scaling functions are universal, have a regular expansion in powers of g(\nL), and are 
normalized such that C R (0) = R\ s , C x (0) = Cd R (0) = 1. The additional corrections due to 
the irrelevant operators decay as powers of 1/L. 

The large-L behavior of f3f(L) follows from Eq. (70). Since k(x) ~ x + 0(x 2 ), we obtain 



c^(lnL) 1 / 2 

Pf - Pc= j 



L\/\n L/Lq 



1 _ 6 3 In In L/L ^ q 



2 lnL/Lo \\nL/L 



(74) 



where L is computed at the critical point t — 0. 

We finally mention that Eqs. (65), (66), (67) hold at fixed R^ = Ri s as well. The 
corresponding universal scaling functions depend on the values of U22 at R^ = R\ s fixed, i.e., 
U 22 (L) = U 22 ((3f(L),L), (we denote them by f R (U 22 ), f x (U 22 ), and f dR (U 22 ), respectively) 
and have a regular expansion in powers of U 22 . 
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TABLE I: MC data at fixed R$ = R Is = 0.9050488292(4). For each model and lattice size L, we 
report the number of samples N s , the quartic cumulants U4 and U22, the magnetic susceptibility 
X, the derivative Rl = dR^/dfi, and the specific heat Ch- If the asymptotic behavior is controlled 
by the Ising fixed point, for L — ► 00 we should have U4 — > ?7i s = 1.167923(5) and U22 — ► 0. 



Model 




L 


AV10 3 


Ua 


U22 


X 




c h 


RSIM, p = 


0.9 


8 


5361 


1.16476(3) 


0.05083(3) 


36.1853(9) 


6.5911(11) 


2.7285(5) 






16 


2560 


1.16463(4) 


0.04170(4) 


122.367(4) 


12.608(5) 


3.4282(12) 






32 


1280 


1.16507(4) 


0.03618(4) 


412.573(14) 


24.029(9) 


4.0283(14) 






64 


640 


1.16563(6) 


0.03237(6) 


1389.57(8) 


45.91(3) 


4.550(3) 






128 


640 


1.16597(6) 


0.02947(6) 


4677.0(3) 


87.84(7) 


5.014(3) 






256 


653 


1.16619(5) 


0.02704(5) 


15741.3(7) 


168.50(9) 


5.431(3) 






512 


633 


1.16656(4) 


0.02522(5) 


52962(2) 


324.18(19) 


5.815(3) 


RSIM, p = 


0.7 


8 


640 


1.14557(10) 


0.09561(10) 


25.841(3) 


1.2536(13) 


0.2155(3) 






16 


2176 


1.15206(6) 


0.07526(6) 


86.941(5) 


2.6583(11) 


0.30474(14) 






32 


1280 


1.15682(6) 


0.06297(7) 


293.888(15) 


4.967(2) 


0.35203(12) 






64 


658 


1.15996(7) 


0.05491(8) 


993.26(6) 


9.140(4) 


0.38283(10) 






128 


843 


1.16185(6) 


0.04871(7) 


3351.92(14) 


16.903(6) 


0.40516(7) 






256 


1288 


1.16313(4) 


0.04368(5) 


11299.1(3) 


31.501(9) 


0.42262(5) 


± J Is, p = 


0.95 


8 


3200 


1.16399(2) 


0.04026(3) 


40.9962(8) 


6.0887(15) 


2.6027(7) 






16 


3200 


1.16405(3) 


0.04023(3) 


139.318(5) 


11.163(2) 


3.1527(8) 






32 


3200 


1.16439(2) 


0.03847(3) 


470.511(11) 


20.924(3) 


3.6299(5) 






64 


812 


1.16482(4) 


0.03592(5) 


1585.27(7) 


39.570(10) 


4.0371(7) 






128 


658 


1.16527(4) 


0.03331(6) 


5337.0(3) 


75.12(2) 


4.3865(5) 
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IV. FINITE-SIZE SCALING ANALYSIS OF MONTE CARLO SIMULATIONS 



A. Monte Carlo simulations 

We perform high-statistics MC simulations of the RSIM at p — 0.9, 0.7, and of the ± J 
Ising model at p — 0.95. We consider square lattices of linear size L with periodic boundary 
conditions. In the MC simulations of the RSIM we use a mixture of Metropolis and Wolff 
cluster 61 updates as we did in the three-dimensional case reported in Ref. 42. In the case of 
the ±J Ising model, the Wolff cluster update is expected to be slow 62 so that we only use 
Metropolis updates with multispin coding. 

Instead of computing the different quantities at fixed Hamiltonian parameters, we com- 
pute them at fixed = £/L = R\ s . This means that, given a MC sample generated at 
P — A-un, we determine the value f3f such that R^((3 = (3f) = Rf. All interesting observables 
are then computed at (3 — /3f. The pseudocritical temperature f3f converges to j3 c as L — > oo. 
This method has the advantage that it does not require a precise knowledge of the critical 
value P c (an estimate is only needed to fix /3 run that should be close to (3 C ). Moreover, for 
some observables the statistical errors at fixed R^ are smaller than those at fixed /3 (close 
to j3 c ). A2,m In order to compute any quantity at ft — f3f, we determine its Taylor expansion 
around /3 run , as we did in our previous work. 62 Particular care has been taken to avoid any 
bias due to the finite number of iterations for each sample: we use the method described in 
Ref. 42 and extended to correlated data in Ref. 62. 

The results at fixed R% = R is are reported in Table I. For each model and lattice size 
L, we report the number N s of samples, the MC estimates of the quartic cumulants U4 
and U22 at fixed R% = R\ s (we denote them with U4 and U22, respectively), the magnetic 
susceptibility Xi the derivative R'^ = dR^/df3, and the specific heat Ch- 

B. Results 

1. Approach to the 2D Ising fixed-point values 

Since we perform our FSS analysis keeping R^ = i? Is fixed, if the critical behavior is 
controlled by the Ising fixed point, in the large-L limit we should have 

U 22 (L) -> 0, U 4 (L) -> U ls , (75) 
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FIG. 2: The phenomenological coupling U22 vs 1/lnL. The lines show the results of fits to 
Eq. (87). For the RSIM at p = 0.9 and the ±J Ising model we fit all data, while for the RSIM 
at p = 0.7, we use data satisfying L > 16. Note that the asymptotic slope as 1/lnL — ► of 
the resulting curves is identical in the three cases, confirming the universality of 6*22,1, defined by 
U22 = C22,i<?(hiL) + 0(5(lnL) 3 ), see Sec. IV B 3. 

where 54 U\ s = 1.167923(5) is the universal large-L limit of the quartic (Binder) cumulant 
at the critical point in the 2D Ising model. Since disorder is expected to be marginally 
irrelevant, see Sec. Ill A, the approach of U22 and U4 to their large-L Ising limit is expected 
to be logarithmic. 

The MC data of U4 and U22, reported in Table I, clearly approach the Ising values (75). In 
the case of U4, see Table I, the MC data are very close to Ui s = 1.167923(5). For the largest 
lattices the relative difference S 4 = \U 4 - U la \/Ui s is very small, 5 4 w 0.0012, 0.0041, 0.0023 
for the RSIM at p = 0.9 (L = 512) and p = 0.7 (L = 256), and the ±J Ising model at 
p = 0.95 (L = 128), respectively. However, the asymptotic approach to the large-L Ising 
value is very slow, hinting at logarithmic corrections. This is also strongly suggested by the 
MC data of U22, which are shown versus 1/lnL in Fig. 2. 

In order to check the approach of the critical exponents to the Ising values, we define the 
effective exponents 



VeS (L) = 2 - 



ln X (2L)/ X (L) 
In 2 



and 



I ML) = 



\nR' $ (2L)/R'^L) 



In 2 



(76) 



(77) 
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FIG. 3: MC estimates of r] e g(L). The dashed line corresponds to the Ising value r\ = 1/4. The 
dotted lines are drawn to guide the eye. 
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FIG. 4: MC estimates of 1/v c k{L). The dashed line corresponds to the Ising value 1/v = 1. The 
dotted lines are drawn to guide the eye. 



l/v UeS (L) 



]nUj(2L)/Uj(L) 
In 2 



(78) 



where we indicate the derivative with respect to f3 with a prime. The MC estimates of 
rjes(L) and l/u e s(L) are plotted in Figs. 3 and 4. They appear to approach the Ising 
values 7] = 1/4 and 1/v = 1. In the case of t], the raw data are already very close to the 
Ising value: the largest lattices give r/ eff (L = 256) = 0.24959(8) for the RSIM at p = 0.9, 
r] cS (L = 128) = 0.24686(8) for the RSIM at p = 0.7, and r? eff (L = 64) = 0.24871(10) for 
the ±J Ising model at p — 0.95. In the case of l/z/ e ff(L) the approach is much slower. At 
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the largest available lattices we find l/V eff (L = 256) = 0.9441(12) for the RSIM at p = 0.9, 
l/uas(L = 128) = 0.8981(7) for the RSIM at p = 0.7, and l/is eS (L = 64) = 0.9249(5) for 
the ±J Ising model at p — 0.95. Anyway, all data show an upward trend towards the Ising 
value \ jv — 1 . 

These results provide already a quite strong evidence that the asymptotic behavior of the 
FSS is universal and it is controlled by the Ising fixed point, with scaling corrections which 
decay very slowly. In the following we report a more careful analysis of these logarithmic 
corrections, showing that they have a universal pattern which is consistent with the RG 
predictions obtained in Sec. III. 



2. Universal finite-size behavior as a function of the phenomenological coupling U22 

As discussed in Sec. Ill B, the FSS formulas obtained from the RG equations of Sec. Ill A 
can be written in terms of the phenomenological coupling U22- In the following we compare 
the MC data with the predictions reported in Sec. Ill B and III C, and in particular with 
Eqs. (65), (66), and (67). 

Let us first consider the quartic cumulant U 4 defined in Eq. (14). At fixed U 4 (L) is 
expected to behave as 

U 4 (L) = f U4 (U 2 2), (79) 
where fu 4 (%) is a universal function, analytic at x — 0, satisfying /[/ 4 (0) = U\ s . Corrections 
to the behavior (79) vanish as powers of 1/L. The scaling behavior (79) is well satisfied by 
the MC data, as shown in Fig. 5. All data fall on a single curve, except for a few of them 
corresponding to small values of L (this is particularly evident in the data for the RSIM at 
p = 0.9), indicating the presence of power-law scaling corrections. The results show that 
the linear term is absent or negligible in the expansion of fu 4 {U22) around U22 — 0; if the 
data are plotted versus C/| 2 , they fall on an approximately straight line, suggesting that 

U 4 (L) - U u = cU 2 2(L) 2 + 0{Ul 2 ). (80) 

A fit of the numerical results to U 4 (L) — Uj s = c£7 22 (£) 2 gives c = 2.4(2). This implies that 

^ L) - u '- + wm^ + - (81) 

where L Q is the model-dependent constant that appears in the expansion of <?(mL) (as such, 
it is independent of the quantity that one is considering). 
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FIG. 5: Uja - U A vs U22 (above) and C/| 2 (below). 
As discussed in Sees. Ill B and III C, at fixed R^, x behaves as 

X = d x L 7 l i f x (U 22 {L)), 



(82) 



where f x (x) is a universal function such that f x {0) = 1. This means that, by properly 
choosing constants e x , the combination e x xL~ 7 / 4 is a universal function of U^. In Fig. 6 we 
show this quantity. The plot is clearly consistent with Eq. (82). Note also that if the data are 
plotted versus [7f 2 they approximately fall on a straight line, suggesting f x (x) = 1 + 0(x 2 ), 
analogously to the case of U4. 

In Fig. 4 we showed the effective exponents (77) and (78) related to the thermal exponent 
v. The data approached the Ising value P\ a = 1 with slowly decaying corrections. The 
effective exponents computed by using Eqs. (77) and (78) were very close, as shown in Fig. 4 
for the RSIM at p — 0.9 (this is also true for the other model considered). For this reason, 
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FIG. 6: Plot of ln(e x xL- 7 / 4 ) vs U 22 (top) and vs C/| 2 (bottom); we set e x = 1,1.4,0.88 for the 
RSIM at p = 0.9 and p = 0.7, and for the ±J Ising model. The constants e x have been chosen 
such as to obtain the best collapse of the MC data. 

in the following we focus on R'^. As discussed in Sec. Ill B and III C, the derivative R'^ at 
fixed R^ scales as 

4 = d dR LU 22 (L) 1 / 2 f dR (U 22 (L)), (83) 

where fdnix) is a universal function. This means that, by properly choosing some constants 
edR, the combination e^i?^/ ' L is a universal function of U 22 . In Fig. 7 we show such quantity. 
The plot is clearly consistent with Eq. (83): the data fall on a single curve and approach 
zero as U 22 {Lfl 2 when U 22 -»• 0. A gain, note the presence of power-law corrections for large 
values of U 22 . 

The approach of f3f(L) to (3 C is given by Eq. (74). Equivalently, we can also consider the 
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FIG. 7: Plot of e dR R'jL versus Tj\£ '. We have chosen e dR = 1, 6.9, 1.2 for the RSIM at p = 0.9 
and p = 0.7, and for the ± J Ising model. The constants e^R have been chosen such as to obtain 
the best collapse of the MC data. 



scaling form 



(84) 



L L 

We determine /3 C by performing fits to Eq. (84). We include only data such that L > L min , 
where L m i n is the smallest cutoff which provides fits with % 2 /DOF < 1. Moreover, as a 
check we have also performed fits to Eq. (84) in which we only consider the leading term 
(i.e. we set d 2 = 0). For the RSIM at p = 0.9 we obtain f3 c = 0.525838(1) (fit with d 2 = 0) 
and j3 c = 0.525835(2), if both terms are taken into account. Analogously, these two fits 
give p c = 0.93294(1), 0.93289(3) for the RSIM at p = 0.7 and (3 C = 0.53362(1), 0.53348(2) 
for the ±J Ising model. Our final estimates are (3 C = 0.525835(3), 0.93289(5), 0.5335(1) 
respectively for the RSIM at p — 0.9 and p = 0.7, and the ±J Ising model at p = 0.95. 
Consistent results are obtained by fitting the data of (3f(L) to Eq. (74). 



3. Universal logarithmic corrections as a function of L 



In the following we directly check the dependence on L of U22, R'^, and of the specific 
heat Ch- As discussed in Sec. Ill C, for L — > 00 the phenomenological coupling U22 behaves 
as 

U2 2 (L) = C 22A g(\nL)+0(g 3 ), (85) 
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where 622,1 is a universal constant. The absence of the term of order g 2 fixes uniquely the 
normalization of the coupling g. This quantity can be expanded in powers of 1/ ln(L/L ) to 
different orders. Using the expansion (34) with y = In L/L , we can perform three different 
types of fit, corresponding to three different approximations for g(lnL). In fit (a), we fit 
U 22 {L) to 

= nrfe (86) 

where 6*22,1 and L are free parameters. In fit (b), we also include the next term proportional 
to 63, i.e. we fit the data to 

u (t\ ^2,1 6 22 ,i& 3 In In L/L 

U 22 {L - — — J T IT \2 ' 

InL/Lo (InL/Lo) 2 
where 622,1, &3, and L are free parameters. Finally, we can also include the next term 
obtaining [fit (c)] 

J 1 6 3 lnlnL/L fe 2 [(lnlnL/L ) 2 - lnlnL/Lp] 1 

W) " C ^ [hTL/Lo - (lnL/L ) 2 + J ' ( } 

where 622,1, b 3 , and L are free parameters. The results of the fits for different values of L min 
are reported in Table II. Let us consider first the fit of the data for the RSIM at p — 0.9 for 
which we have the largest lattices. Fit (a) has a very poor x 2 , indicating that the data are 
not well fitted by a single logarithmic term. If we include the next correction the x 2 drops 
dramatically, indicating that our results are precise enough to be sensitive to the elusive 
lnlnL/L terms. Beside the very good x 2 , the results are also very stable with L min . This 
stability should not be trusted too much however, since fit (c) — which a priori should better 
since we include an additional set of corrections — has a very poor x 2 and gives results that 
vary somewhat with L min . As an additional check we also fit IJ^iiJ-) to 

U 22 (L) = C 22 yg(\nL) + C 22 , 3 g(\nL) 3 , (89) 

using the expansion of g(lnL) used in fit (c). For L mhl = 8,16,32 we obtain x 2 /DOF 
= 213/3, 135/2, 16/1; they are better than those obtained in fit (c), but still sig- 
nificantly worse than those obtained in fit (b). Correspondingly, we obtain 622,1 = 
0.210(1), 0.254(3), 0.268(10) and 6 3 = 1.44(1), 0.89(10), 1.0(3) for the same values of L min . 
Finally, we fit U 22 (L) to Eq. (89) by using the exact expression for ^(lnL): for each L/L , 
g(\nL) is obtained by inverting F(g) = In L/L , where F(x) is defined in Eq. (32). If we 
only include the leading term, i.e. we set 622,3 = 0, the quality of the fit is significantly 
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worse than that of fit (b) and better than that of fit (c): x7 D0F = 515/4,52/3,6/2 for 
L min = 8, 16, 32. Correspondingly, C 22 ,i ~ 0.27, 0.29, 0.31 and b 3 « 2.0, 2.4, 2.7. Though the 
scatter of the estimates of C 2 2,i is significantly larger than the statistical errors — this should 
be expected since x 2 /DOF is significantly larger than 1 in most of the cases — the data show 
a clear pattern. If we take the central estimate from fit (b), we obtain C 2 2,i ~ 0.28. To 
estimate a reliable error, note that all results of the fits with L m i n > 16 lie in the interval 
0.23 < C22,i < 0.31. A conservative error is therefore ±0.05, so that C*22,i — 0.28(5). It 
is more difficult to estimate 63, since this parameter varies significantly from one fit to the 
other. In any case, note that all results satisfy 63 > 0, in contrast with the theoretical 
prediction 63 = —1/2. It is not clear if this difference should be taken seriously. It might 
be that it is only an indication that we are not sufficiently asymptotic to estimate correctly 
the coefficient of the slowly varying lnlnL/Lo term. 

Since C22,i is universal, we can check its estimate by comparing the above-reported results 
with those obtained in the other two models, for which we have less data. For both models, 
fit (a) is significantly worse than fit (b) or fit (c). For the RSIM at p = 0.7 fit (b) and fit 
(c) have similar reliability. The corresponding estimates of C 2 2,i are fully consistent with 
that reported above. For the ±J Ising model, only fit (c) is reliable. The estimates of 
6*22,1 are again consistent with those obtained in the RSIM. The universality of the leading 
logarithmic correction is well satisfied by our data. 

The scale L is very poorly determined and varies significantly with L min and the type of 
fit. The ratio of the scales can also be determined by directly matching the numerical data. 
If power-law scaling corrections are negligible, we should have 

k22,modcl l{L) = ^ r 22,modcl2( K -^) (90) 

for some constant k, which is the ratio of the scales L pertaining the two models. By direct 
comparison we obtain L ,rsim, p =o.7 ~ ^o,rsim, p =o.9, k > 16, and L ,±j ~ «^o,rsim, p =o.9, 
with 2 < k < 4. Since L is independent of the observable, these relations should not be 
specific of U22 but should apply to any RG invariant quantity: indeed, as can be seen from 
the data reported in Table I, they also approximately hold for U4. Note that Lq increases 
with p in the RSIM as expected: the Ising critical behavior is observed for L > L min , with 
L min increasing with p. 

In order to check the L-dependence of the derivative R'^, previosly discussed as a function 
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TABLE II: Results of the fits. We do not report the results of fit (b) for the ± J Ising 
£min = 16 because this fit is unstable (apparently, the \ 2 continuously decreases as 63 
Lq — > 0). DOF is the number of degrees of freedom of each fit. 



model with 

— ► — 00 and 





Fit (a) 




Fit (b) 






Fit (c) 




-^min 


y 2 /DOF 

A. / 


C22 1 


X 2 /DOF 


C22,l 


b 3 


X 2 /DOF 


C*22,l 


b 3 










RSIM p 


= 0.9 








8 


1844/5 


0.193(1) 


3.4/4 


0.280(2) 


1.35(1) 


1280/4 


0.222(1) 


0.91(3) 


16 


221/4 


0.227(1) 


3.1/3 


0.281(3) 


1.36(3) 


164/3 


0.240(2) 


0.85(7) 


32 


27/3 


0.235(2) 


3.1/2 


0.281(5) 


1.36(8) 


20/2 


0.252(5) 


0.88(23) 










RSIM p 


= 0.7 








8 


748/4 


0.276(1) 


15/3 


0.356(2) 


0.88(2) 


37/3 


0.334(1) 


1.09(1) 


16 


95/3 


0.287(1) 


14/2 


0.350(5) 


0.83(5) 


1.7/2 


0.324(1) 


1.30(1) 


32 


0.6/2 


0.297(1) 


0.3/1 


0.28(3) 


-0.3(7) 


0.3/1 


0.28(3) 


-0.3(5) 










±J model 








8 


4211/3 


0.986(4) 


2753/2 


0.610(4) 


2.01(1) 


27/2 


0.345(3) 


1.90(2) 


16 


389/2 


0.449(4) 








0.02/1 


0.315(6) 


1.79(2) 



of £7 2 2, we perform fits of the MC data of R'^ to the behavior 



, R'e , , L a 3 lnlnL/L a 4 
m ~T ~ °l m m T \~ Q>2 -\ : — h ; — , 



(91) 



L Lq \tlL/Lq \hL/Lq 

taking; and Lq as free parameters. According to theory, we should find ai = —1/2. 

Because of the presence of 5 free parameters this fitting form can be safely used only for 
the RSIM at p = 0.9. If we fit all data we obtain a t = -0.53(5) (x 2 /DOF = 0.71); if we 
discard the result corresponding to L = 8 we obtain ai = —0.44(11). These results are in 
good agreement with the theoretical prediction a\ = —1/2. 
Finally, we consider the specific heat Ch, 

\{n 2 ) - (h) 2] 



V 



(92) 



where H is the Hamiltonian. The RG analyses of Refs. 1,2,4,28 predict a diverging lnlnL 
asymptotic behavior. In Fig. 8 we show the MC data of Ch at fixed R^ = Ri s . They are 
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a RSIM p=0.7 
_ ■ +-J p=0.95 
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0.8 1.0 1.2 1.4 1.6 1.8 

In InL 

FIG. 8: MC data of the specific heat versus In InL. The lines correspond to fits to Eq. (93). 

definitely consistent with the theoretical prediction for its asymptotic behavior 

C h « Alnln(L/L Q ) + B, (93) 

where A and B are nonuniversal parameters, and L is the model-dependent scale. Fits of 
the data to Eq. (93), taking A, B, and L as free parameters, suggest A ~ 5, 0.1, 2 for the 
RSIM at p = 0.9 and p = 0.7, and the ±J Ising model at p = 0.95, respectively. There is a 
large difference between the values for RSIM at p = 0.9 and p = 0.7, but this should not be 
surprising because the p = 0.7 is quite close to the percolation point 46 p = p perc ~ 0.59. 
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